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ABSTRACT 

The strong dependence of the neutrino annihilation mechanism on the mass accretion 
rate makes it difficult to explain the LGRBs with duration in excess of 100 seconds 
as well as the precursors separated from the main gamma-ray pulse by few hundreds 
of seconds. Even more difficult is to explain the Swift observations of the shallow 
decay phase and X-ray flares, if they indeed indicate activity of the central engine 
for as long as 10 4 seconds. These data suggest that some other, most likely magnetic 
mechanisms have to be considered. Since the efficiency of magnetic mechanisms does 
not depend that much on the mass accretion rate, the magnetic models do not require 
the development of accretion disk within the first few seconds of the stellar collapse 
and hence do not require very rapidly rotating stellar cores at the pre-supernova state. 
This widens the range of potential LGRB progenitors. In this paper, we re-examine 
the close binary scenario allowing for the possibility of late development of accretion 
disks in the collapsar model and investigate the available range of mass accretion rates, 
black hole masses, and spins. We find that the black hole mass can be much higher 
than 2 — 3M Q , usually assumed in the collapsar model, and normally exceeds half of 
the pre-supernova mass. The black hole spin is rather moderate, a = 0.4 — 0.8, but 
still high enough for the Blandford-Znajek mechanism to remain efficient provided the 
magnetic field is sufficiently strong. Our numerical simulations confirm the possibility 
of magnetically driven stellar explosions, in agreement with previous studies, but point 
towards the required magnetic flux on the black hole horizon in excess of 10 28 Gcm 2 . 
At present, we cannot answer with certainty whether such a strong magnetic field can 
be generated in the stellar interior. Perhaps, the supernova explosions associated with 
LGRBs are still neutrino-driven and their gamma-ray signature is the precursors. The 
supernova blast clears up escape channels for the magnetically driven GRB jets, which 
may produce the main pulse. In this scenario, the requirements on the magnetic field 
strength can be lowered. A particularly interesting version of the binary progenitor 
involves merger of a WR star with an ultra-compact companion, neutron star or black 
hole. In this case we expect the formation of very long-lived accretion disks, that may 
explain the phase of shallow decay and X-ray flares observed by Swift. Similarly long- 
lived magnetic central engines are expected in the current single star models of LGRB 
progenitors due to their assumed exceptionally fast rotation. 

Key words: black hole physics - accretion disks - supernovae: general - gamma-rays: 
bursts - binaries: close - MHD - relativity 



1 INTRODUCTION 

The nature of Gamma Ray Bursts (GRBs) remains one of 
the most intriguing problems of modern astrophysics. It is 
now widely accepted that the gamma ray emission is gener- 
ated in ultrarelativistic jets but many basic questions related 
both to the physics of these jets and to the mechanisms 
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of their production remain open. Although many promis- 
ing theories have been developed over the years since the 
discovery of GRBs, we are still some way out from solid 
understanding of this phenomenon. For example, the ob- 
served connection between the Long-duration Gamma Ray 
Bursts (LGRBs) and supernovae (SNe) indicates that these 
bursts are connected to deaths of massive stars but the 
details are not clear. In one model of LGRBs, the stellar 
collapse results in a normal successful supernova explosion 
but the newly born neutron star (NS) is very unusual. It 
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has both exceptionally high magnetic field, and for this 
reason it is called a magnetar, and extremely rapid rota- 
tion f e.g.lUsovlll992l;l Thompson et al. ll2004l : lMetzger et al. I 
120071 ; lUzdenskv fc MacFadvenll2007l ). The powerful magne- 
tohydrodynamic (MHD) wind produced by such remnant is 
capable of both accelerating the supernova shell above the 
expansion speed of normal supernovae, to the level of hyper- 
nova e, and production of collimated ultra-relativistic pol ar 
jets (iKomissarov fc Barkov I [20071 ; IrBucciantini et al.l l2009). 

In another model, the normal supernova explosion fails 
and the proto-neutron star promptly collapses into a black 
hole (BH). However, the rapid rotation of the stellar pro- 
genitor prevents the rest of the star from falling directly 
into the black hole and a massive neutrino-cooled accre- 
tion disk is formed instead. This allows to turn the failed 
supernova into a successful stellar explosi on, as this disk 
can release enormous amoun ts of energy (IWooslev 1 1 19931 : 
iMacFadven fc Wooslevl [l999) . One way of "utilising" this 
energy is via the neutrino- or magnetically-driven wind 
from the disk. Such wind is not expected to be relativis- 
tic due to the high mass loading at its base. However, 
the polar region just above the black hole is less likely to 
become mass-loaded by the disk matter and can become 
relativistically hot via annihilation of neutrinos and anti- 
neutrinos emitted by the disk. This opens a possibility of 
drivi ng ultra-relativistic LGRB jet s in the collapsar s cenario 
(e.g. IMacFadven fc Wooslevl fl999l : lAlov et aljlioool l. How- 
ever, the efficiency of this type of neutrino heating is a very 
strong function of both the mass accre tion rate and the ro- 
tation rate of the central black hole ([ Popham et al. 1 19991 : 
IChen fc Beloborodovll2007l : IZalamea fc Beloborodovll 2009h . 



According to the calculations by iPopham et al 



19991 ) for 

the black hole with the rotation parameter a — 0.5 the 
energy deposition rate via the neutrino annihilation pro- 
cess drops from L vi ? = 4 x 10 48 erg/s for M = 0.1M Q to 
L V v = 6 x 10 44 erg/s for M — 0.01M Q and for the accretion 
rate of M = 0.1M Q from L vB = 2 x 10 51 erg/s for a = 0.95 
to L V y = 3 x 10 48 erg/s for a = 0. Therefore, this version 
of the collapsar model, similarly to the magnetar model, re- 
quires very rapid rotation of the stellar core prior to the 
collapse so that the accretion disk is formed early on, when 
the accretion rate is still high e nough, and t he BH is born 
rapi dly rotating. T h e resu lts by iBirkl et all (|2007l ) suggest 
that IPopham et "all (| 19991 ) may have overestimated the effi- 
ciency of neutrino mechanism for high a. This is because the 
energy released by the disk powers not only the outflow but 
also the flow into the BH. As a increases the inner bound- 
ary of the disk moves closer to the BH and a larger fraction 
of the total neutrino-antineutrino annihilation occurs in the 
region where the vector of dep osited mom e ntum points to- 
wards the black hole. In fact, IBirkl et all (120071 ) find that 
the efficiency of the neutrino annihilation mechanism peaks 
at a ~ 0.6. 

It turns out that such a fast rotation cannot be a general 
result of stellar evolution. Although young massive star often 
rotate sufficiently rapidly at birth, their cores are expected 
to experience strong spin-down during the red giant phase 
and during the intensive mass loss peri od characteristic fo r 
massive stars at the Wolf-Rayet phase (|Heger et al. | [2005). 
In fact, this theoretical result agrees very well with the ob- 
served rotation rates of newly born pulsars. Thus in order 
to retain the rotation rate required in the collapsar model, 



the evolution of LGRB progenitors must proceed along a 
rather exotic route. Recently, it was proposed that a com- 
bination of low metalicity and extremely fast initial rota- 
tion, at arou nd 50% of the break-up speed, could lead to 
such a route |Yoon fc Langerll2005l : IWooslev fc Heger!l2006l : 
lYoon et al . 2006). On one hand, the mass-loss rate decreases 
significantly with metalicity, leading to a significant reduc- 
tion in the total loss of angular momentum. On the other 
hand, the rotationally induced circulation becomes very ef- 
fective at such a high rotation rate and may result in chem- 
ically homogeneous stars that avoid the development of ex- 
tended envelops and hence the spin-down of stellar cores 
via interaction with these envelopes. Moreover, the star re- 
mains compact by the time of its collapse so the LGRB jet 
can break out from the star on the time scale compatible 
with the observed durations of LGRBs. 

Another exotic scenario involves close high-mass bi- 
nary systems, where the fast rotation of stellar cores is 
sustained via the tidal interact i on b e tween companions 
jTutukov fc Cherepashchuk||2003l. |2004 llzzart et al 1 12004 



iPodsiadlowski et alj|2004 Ivan den Heuvel and YoorfeoTl ) 
In this case, the pre-supernova is a compact helium star, 
essentially a Wolf-Rayet (WR) star, because the extended 
envelope is dispersed into the surrounding space during the 
common envelope phase. The stellar rotation in such systems 
is synchro nised with the orbital motion on a very short time 
scale fe.g. lvan den Heuvel and Yoonll2007l ). The contraction 
of CO-cores during stellar evolution leads to their additional 
spin-up but due to the core-envelope coupling only a frac- 
tion of their angular momentum is retained (|Yoon et al.l 
2006). As the result, the core rotation rate is insufficient 
in the cases where the compani on of the helium star is a 
main- sequence star. According to Ivan den Heuvel and Yoonl 
|2007l ). the core rotation can be high enough to fit the col- 
lapsar model with the neutrino-driven LGRB jet only if 
the component is also a compact star, namely NS or BH. 
Three examples of such systems are known to date: Cyg 
X-3, IC 10 X-l, and NGC 300 X-l. Cyg X-3 has a ver y 
short orbital period, only 4.8h l|van Kerkwiik at~aT1ll992l ). 
and the r adius of the WR star in this system is less than 
3 - 6i? jCherepashchuk fc Mofatlll994T ). The recently dis- 
covered IC 10 X-l and NGC 300 X-l have the orbita l 
periods of 35h and 33h respectively dCarpano et"all 120071 : 
IPrestwich et~ell 120071 : ISilverman fc Filippenkol 120081 ). The 
masses of Wolf-Rayet stars are estimate d at 18 — 40Mp) for 
NGC 300 X-l and ~ 35M for IC 10 X-l jClark fc Crowtherl 
|2004|). Given the observed producti on rate of such sys- 
tems |ran~^«njIeuveriind~^^on| l|2007l ) predicted one hyper- 
nova/LGRB every 2000 years in a galaxy similar to our own. 

The neutrino heating is not the only possible 
mechanism behind the explosions of collapsing stars. 
Perhaps somewhat less popular, but the magnetic 
mechanisms are also regard ed as potentially impor- 



tant ( e.g. Bisnovatvi-Koganl Il970l: LcBlanc & Wilson I 



197d: IMacFadven et al.l 120011 : iMoiseenko et al. ] 120061 : 



Burrows et al. 1 120071 ). Likewise, the LGRB jets can also 
be powered via a magnetic mechanism, in particular 
the Blandford-Znajek mechanism, which utilises the ro- 
tatio n al energy of the BH (e.g. [B landford fc Znaiek 



1 19771 : iMeszaros fc Reed Il99?l: iLee et all l2000l: iProga et al 
20031: iMcKinnevI 120061: iBarkov fc Komissarov I l2008al : 



Komissarov fc Barkov 2009). Here the black hole is also re- 
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quired to rotate quite rapidly. However, the efficiency of this 
mechanism is not that sensitive to the mass accretion rate 
and such rapid rotation does not have to be achieved right 
after the collapse of the Fe core. Instead, it can be built 
up gradually during the rest of the stellar collapse. This 
difference in the sensitivity to mass accretion rate favours 
the BZ mechanism over the neutrino mechanism in the case 
of very long-duration L GRBs, more than 100 seconds long 
(|MacFadven et al.ll200lT ). The discovery by Swift of the shal- 
low deca y phase and l ate flares in the X-ray light curves of 
LGRBs jzhandl2007l ; IChincarini et all 120071 ) also suggests 
that the central engine may remain active for as long as 
10 4 seconds (e.g. iLipunov fe Gorbovskovll2007l . l200Sf ). Since 
the neutrino mechanism requires the mass accretion rate to 
stay above fewxlO _2 M0/s, such a prolonged activity im- 
plies the progenitor mass in excess of fewxlO 2 M0, which is 
highly unlikelj0- 

Another problem for the model of neutrino-driven 
GRBs are the strong precursors sometimes observed be- 
fore t he arrival of the main gamm a-ray pulse dBurlon et al. I 
120081 ). According to the analysis of lWang fc Meszarosl (120071 ) 
such precursor and the main pulse can be attributed to a 
single eruptive event only when the precursor and the main 
pulse are separated by few seconds. However, in some GRBs 
the delay can be as long as few hundreds of seconds and in 
such cases it is much more likely that the precursor and the 
main pulse correspond to two different events in the life of 
the central engine. They proposed that the precursor is pro- 
duced during the supernova explosion, in the jet powered by 
a rotating magnetised neutron star, and that the main pulse 
is produced during the fallback phase when the neutron star 
collapses into a black hokfl The typical mass accretion rates 
in the fallback scenario, lO^-lO^Mos -1 , are too low for 
the neutrino annihilation mechanism and thus this expla- 
nation implies magnetic origin for the main pulse as well 
|MacFadven et al.ll200ll ). 

Thus, the observations require to include the magnetic 
mechanism, either in the black hole or, in fact, in the disk 
version, or both, in the collapsar scenario. This widens the 
range of potential progenitors of LGRBs. Indeed, we no 
longer need to constrain ourself to the stars with extremely 
rapidly rotating cores but can also include the cases with 
slower rotation where the accretion disk forms much later 
during the course of stellar collapse. 

In this paper we re-examine the scenario of binary pro- 
genitor of LGRBs allowing for the late formation of accretion 
disks and lower mass accretion rates compared to those re- 
quired in the collapsar model with the neutrino mechanism. 
In Section [2] we determine the parameters of binary systems 
which allow formation of accretion disks during the collapse 
of WR companion. We also estimate masses and spins of the 
black holes by the time of accretion disk formation using sim- 
plified analytical model for the structure of pre-supernovae 



1 Typically, the mass of WR star is 9-25 Mq, though some 
observations suggested that it can be as high as 83 Mq 
jSchweickhardt et al.lll999l ; ICrowther 1120071) . 

2 This model may struggle to explain delays shorter than the 
typical fallb ack time, 100-1000 secon ds, found in one-dimensional 
simulations MacFadycn et al. | d200lh . However, due to the rota- 
tional effects the supernova explosions could be highly aspherical, 
resulting in shorter fallback time scales in the equatorial region. 



due to iBethei |l990l ). In Section [3] we investigate the degree 
to which the black hole spin can increase later on, during 
the disk accretio n phase, usin g the s ame approach as in the 
recent study by Ijaniuk et al.l l|2008l ). Here we consider not 
only the Bethe's model but also the polytropic model and 
the models of pre-supernovae based on detailed calculations 
of stellar evolution. In Section [4] we describe the numerical 
simulations of LGRB jet formation with setup based on the 
results obtained in the previous Sections. In Section [S] we 
analyse the potential of the binary scenario in the extreme 
case, which involves merger of the WR star with its ultra- 
compact companion, BH or NS. In Section [S] we summarise 
our main results and discuss their astrophysical implications. 



2 FORMATION OF ACCRETION DISC 

In a synchronised binary the tidal torques force the compo- 
nents to spin with the same rate as the orbital rotation, 



fi 2 =GM s (l + 9 )/L 3 



(1) 



where L is the orbital separation, M s is the mass of the star 
under consideration, and q = M com /M s , where M com is the 
mass of the companion star. Since, the orbital frequency de- 
creases with L, the maximum possible spin is reached when 
the separation is minimum. This corresponds to the case 
where the star radius is about the size of its Roche lobe. 
The relation between the minimum separation L m i n , the 
stellar radius R a , and q can be approximated with sufficient 
accuracy for 1/100 < q < 100 as 



2Mq 02084 R s 



(2) 

l|Plavec fe Krotochvillll964l ). 

During the stellar collapse the centrifugal force will halt 
the free-fall of the outer layers and promote the development 
of accretion disk provided the specific angular momentum 
on the stellar equator exceeds that of the marginally bound 
circular orbit for the black hole with the same mass and 
angular momentum as the star. The angular momentum of 
Kerr black holes is 



Jh 



GMi 



where — 1 < a < 1 is the dimensionless spin parameter and 
Mh is the hole mass. The specific angular momentum of test 
massive particles on circular orbits in the equatorial plane 

is 



(r 2 - 2ar 1 / 2 + a 2 ) 



GM h 



(3) 



r 3/4( r 3/2„ 3r l/2 + 2a )l/2 c ' 

where r = R/R g and R g — GMh/c 2 , and the radius of the 
marginally bound orbit is 

r mb = {2-a + 2(l-a) 1/2 } (4) 
|Bardeen et al.lll972l ). The disk formation condition is 



Q s Rt > I 

n 



(3) 



where l m b = l(r m b)- As we shell see later, at the time of the 
disk formation a is quite small. Using the Taylor expansion 
we find that 



l mb = (4-a)^+0(a 2 ) 



4GMft 



(6) 
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Using this result and Eq[T]we can now write the disk forma- 
tion condition as 



. 1 + 9 
.Rs) 16 rs 



(7) 



where r s — R s /R gs and R gs = GM B /c 2 . For the typical 
parameters of WR stars this amounts to 

1/3 / M s " " 1/3 



L< UR s (l + q y« ( ^ 



V 10M, 



(8) 



The comparison of this result with Eq[2] shows that collapse 
of WR stars in close binaries can indeed lead to forma- 
tion of accretion discs (|Bisnovatvi-Kogan fc Tutukov|[20o3 ; 
iPodsiadlowski et alj|2004h . We can rewrite the above condi- 
tion in terms of the binary period, Tj,, as 

-i / ,-, \ 2 



T b < -T k r] 



-1/2 



48hr 



/ M 3 



\ 10M G 



Re 



(9) 



where Tk is the Keplerian period at R = R s . This up- 
per limit is about five tim es higher than that obtained in 
IPodsiadlowski et alj (|2004T ) who required the disk to form 
immediately after the collapse of iron core. 

Effectively cooled accretion disks remain geometrically 
thin and their inner radius is given by the radius of the last 
stable circular orbit 



(10) 



(11) 



r ms = {3 + Z 2 - [(3 - Z 1 )(3 + Z 1 + 2Z 2 )] 1/2 }, 
where r ms = R ms / R g and 

Z x = 1 + (1 - a 2 ) 1/3 [(l + «) 1/3 + (1 - a) 1 / 3 ], 
Z 2 = (3a 2 + Zl) 1 ' 2 . 

(Bard een et al.l fl9T2h ■ The corresponding specific angular 
momentum, l ms , determines the evolution of the black hole 
spin via the disk accretiorjf]. 

The outer radius of the disk, Rd, is determined by 
the specific angular momentum on the stellar equatoiQ 
l B — fl s R 2 . Assuming that ra = Rd/R gs 3> 1 the angular 
momentum at the outer edge of the disk is simply 

l d = (GM s R d y 1/2 . 

Matching Id and l s , and using Eq[TJ we find that 

r< = r.(l + *)(£)" 
Thus, 

M, 



(12) 



47 



R s 

Rq 



1OM 



(13) 



where L = (1 + g) _1 ' 3 (L/7? s ). For the widest orbital sepa- 
ration which still allows disk formation (see Eq[8l this equa- 
tion gives r d < 17 whereas for the closest one (see Eq[2} we 
have rd < 5 x 10 3 . Thus, the model predicts a wide range 
of accretion disk sizes. Compact disks and the inner regions 
of large disks will cool via the neutrino emission, whereas 

3 For simplicity, we ignore the effects of magnetic torques on the 
evolution of black hole spin. 

4 In fact, various torques operating in the accretion disk change 
the angular momentum and hence t he location of the outer edge 
but this effect is relativel y minor jShakura fc Sunvaev 1 1 1973t 
IChen fc Beloborodovll2007f) . 





Figure 1. The black hole spin (top panel) and mass (middle 
panel) at the disk formation time as functions of the progenitor 
spin, a s , for M c /M s = 1/3 (dash-dotted line), 1/5 (dotted line), 
1/9 (dashed line), and 1/31 (solid line). The bottom panel shows 
the time of disk formation as a function of a s for the same models. 



the outer regions of large disks will remain adiabatic. The 
accretion time of neutrino cooled disks 



t d 



M h 



\ lOAf© 



6/5 



(14) 



|Popham et all 1 19991 ) is significantly less than the free-fall 
time scale 
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Figure 2. Evolution of the black hole mass and spin in the Bethe's model with M c /M g = 1/9. Left panel: The total mass of black hole 
(solid line) and the mass accumulated via the accretion disk (dash-dotted lines) for the progenitor spin a B = 0.33,0.58,1.0,1.7, and 3. 
The higher value of a 3 correspond to the higher fraction of mass processed via the disk. Right panel: The spin parameter a of the black 
hole for the progenitor spin a s = 0.33 (dashed line), 0.58 (dotted line), 1.0 (thick solid line), 1.7 (dot-dashed line), and 3.0 (thin solid 
line). The evolution time is given in the units of the free-fall time (see Eq ll5D . 



t ff « 240 



R 

Rq 



3/2 



M 3 



\ 10M G 



-1/2 



(15) 



The accretion time of large disks can be estimated using the 
a-model for slim disks (5 — Hd/Rd — 0.3) 



t d 



(16) 



-' 1 0.01 J Kw 3 ) ^IOMq 

|Shakura fc Sunvaev 1 179731 ). Thus, with the exception of 
largest disks, the time scale of disk accretion is shorter com- 
pared to the free-fall time scale, and hence the growth rate 
of the black hole mass is given directly by the rate of the 
collapse. 

In order to estimate the mass and rotation rate of the 
black hole at the time of the disk formation one needs to 
know the mass distribution of progenitor at the o nset of 
collap se. Here we adopt the power law model used bv lBethel 
i| 19901 ) in his analytical models of core-collapse supernovae, 



p{R) = pc 



R> R c 



(17) 



where R c is the radius of iron core. Simple integration allows 
us to find the following equations for the mass 

M(R) = A?K Pc Rl \a(R/Rc) (18) 

and the moment of inertia 

1 M(R)R 2 

i{r) -3Hr7r7) (19) 

of the shell between the iron core and the radius R. 

By analogy with the black hole theory it is convenient 
to describe the rotation rate of collapsing star using the spin 
parameter 

- = mi- (2 °) 

In Bethe's model it relates to f2 s via 



Q, — a. 



3GM 3 (1 + t 1 ) 2 
R 2 s c 



In y 3 



(21) 



where y 3 — R 3 /R c , rj = M c /M s , and we ignore the small 
contribution of compact iron core to the total spin of the 
star. The condition (f2]| with q = 1 implies that 



91ny a 



1/2 



5.2 



R 

Rq 



1/2 



10M, 



-1/2 



(22) 



where we used y s — 100. This seems to suggest that the 
stellar collapse may lead to formation of rapidly rotating 
black holes. 

Suppose that the disk is first formed at time t* and that 
by this time the black hole has swallowed the star up to the 
initial radius R — R,. Assuming that the black hole spin at 
this point is low, a, < 1, we have 



n.ig = (4-a.) G ^ + ^ 



(23) 



where M* = M(R„). Using Eqs.JTSJ and J2TJ) this condition 
can be written as the following algebraic equation for y, = 
R* J R c 



2 

V* = 



(4-a»)j/ s 2 
3(1 + rj) 2 a 3 lny 3 



In y* 
lny 3 



+ V 



where 



1 + 3(lny, + T)lny s ) ' 



(24) 



(25) 



This equation is solved numerically and the results are pre- 
sented in FigUl One can see that the disk is formed relatively 
late, with the typical time t* > O.lt//, when more than a 
half of the star has already collapsed into the black hole. 
However, the black hole spin at this moment is relatively 
low, 0.2 < a, < 0.4. 
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Figure 3. The final value of the black hole spin in Bethe's 
model with as a function of the progenitor spin for models with 
M c /M star = 1/3 (dot-dashed line), 1/5 (dotted line), 1/9 (dashed 
line), and 1/31 (solid line). 



3 GROWTH OF BLACK HOLES 

In order to explore the evolution of the black hole spin for 
more sophisticated models of LGRB progenitors, as well as 
its evolution in Bethe's model after the disk formation, one 
can integrate the following system of dynamic equations 



dMh A E-2 (T>\ 

— = ^R p(R), 



the solid body rotatiorjf] and assumed that initially the 
black hole is rapidly rotating, with a = 0.85. They also 
used the model of geometrically thick and radiatively in- 
efficient disk, with the inner edge located at the radius of 
the marginally bound orbit, whereas we use the thin disk ap- 
proximation, which is more suitable for the neutrino-cooled 
collapsar disks. 



3.1 Bethe's model 

Figure [5] shows the typical evolution of the black hole mass 
and spin, as described by Eqs l26l and 1271 for the Bethe's 
model. One can see that the black hole spin increases sig- 
nificantly above the values attained by the time of disk 
formation. Eventually, it reaches the relatively high val- 
ues of a = 0.3 — 0.8, the final spin depending mainly on 
the progenitor spin and less so on the mass fraction of the 
iron core (see Fig[3}. These higher values of a imply higher 
potential efficiency of both the neutrino annihilation and 
the B landford-Znajek mec h anisms of the LGRB jet produc- 
tion ilPopham et all Il999l; IZalamea fc Beloborodovl I 20091 ; 



iBarkov fc Komissarov 112008a 

The total mass accretion rate can be easily derived from 
the mass distribution and the free fall time (see Eq |30jl : 



M — -- 



3 lny s 



1.45 



10Mp 



M©s~ 



(31) 



(26) 



where t is the time since the start of the collapse. As one can 
see in Figure [2j soon after the disk formation the mass accre- 
tion rate becomes dominated by disk. Initially, the rate can 
be rather high but at around t ~ 100s it becomes insufficient 
for the neutrino annihilation mechanism to operate. 



tt/2 

^ = 4nR 2 p(R) j l(R, 



) sin 9d9. 



(27) 



Here, p(R) is the stellar mass density prior to the collapse 
and l(R,8) is the specific angular momentum retained by 
the fluid element, initially located at the point with the co- 
ordinates {R, 8}, by the time it crosses the event horizon. 
This quantity is given by 



/ = 



l(R,6) 

lms(M h , J h ) 



if 
if 



I < l mb (M h ,Jh) 

I > lmb(M h ,J h ) 



(28) 



where l(R,9) is the distribution of the progenitor's angu- 
lar momentum. The initial conditions for Eqs. (|26l27fl corre- 
spond to the iron core of the WR star 



M h {R c ) = Mo, MRc) = 0, 



(29) 



3.2 Stellar evolution models 

Although the Bethe's model provides a reasonable zero- 
order approximation for the structure of pre-supernova 
stars, the more sophisticated models based on numerical in- 
tegration of the equations of stellar evolution yield some- 
what different stellar structure with wealth of finer de- 
tails. Our next results are based on the pre-collapse struc- 
ture of massive zero age main sequence (ZAM S) stars with 
masse s M 8 = 20Mq and 35Mq described in Irleger et al. I 
(2004). Assuming that stars of close binaries lose their ex- 
tended envelopes we cut of the mass distributions beyond 
the C/O core. This results in the progenitors with masses 
M 3 = 6.15M© ( model A ) and M s = 12.88M© ( model B ) 
respectively, and radius R s — 0.3-Rq.Q The moments of in- 
ertia of models A and B are / ~ 0.065M 3 R 2 3 and 0.074M s i?? 
respectively. Given these parameters, Eqs[T]and[2]imply the 



where M c and R c are respectively the mass and the radius 
of the core. When the accretion rate is determined by the 
free-fall time, R and t can be related via 



2R A 



" ~ 9GM(R) ■ (30) 

The same approach has been used in Ijaniuk et al.l 

(2008) in their search for the laws of rotation that would 
fit the collapsar model of LGRBs. They did not consider 



5 Th e solid body rotation law was studied in Ijaniuk fc Progal 
(2008) but it was assumed there that the black hole was non- 
rotating. 

6 This radius is rather small, twice as small compared to 
the m odels of WR stars constructed in ISchaerer fc Maederl 

19921) and 10-20 time smaller compared to th e observed radii 



Cherepashchuk fc Mo fat 1994; Crowther 1120071) , We can offer no 



clear explanation for this discrepancy. Perhaps, the artificial "re- 
moving" of extended H/He envelope is not a particularly accurate 
procedure. 
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Figure 4. Evolution of black hole's mass and spin in numerical model B. Left panel: The total mass of black hole (solid line) and 
the mass accumulated via the accretion disk (dash-dotted lines) for the progenitor spin a s = 0.20,0.33,0.57,0.96, and 1.6. The higher 
value of a s correspond to the higher fraction of mass processed via the disk. Right panel: The spin parameter a of the black hole for the 
progenitor spin a a = 0.20 (dashed line), 0.33 (dotted line), 0.57 (thick solid line), 0.96 (dot-dashed line), and 1.6 (thin solid line) 
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log 10 t log 10 t 

Figure 5. The accretion rate (in the units of Mq s _1 ) for model A (left panel) and model B (right panel). The solid lines show the 
total accretion rate whereas the dash-dotted lines show the disk accretion rates for different spins of the progenitor; a a = 0.32, 0.54, 0.90, 
1.52, and 2.55 for model A and a a = 0.20, 0.33, 0.57, 0.96, and 1.6 for model B. Higher values of a s correspond to earlier formation of 
accretion disk and higher disk accretion rates. 



spin parameters a s < 2.6 and a s < 1.7 for the models A 
and B respectively; somewhat smaller than in the Bethe's 
model. 

Figure [4] shows the evolution of the black hole's mass 
and spin in model B for different assumed values of the pro- 
genitor's spin. The comparison with the results obtained for 
Bethe's model shows only relatively minor differences, sug- 
gesting that Bethe's model is quite accurate. Figure [5] shows 
the accretions rates, both for the disk and in total, for dif- 
ferent progenitor spins in models A and B. One can see that 
initially the disk accretion rate grows rapidly and soon it 
accounts for most of the total accretion rate. Then it begins 
to decay, approximately as t^ 1 in model A and t~ s in Model 
B. For the cases with faster stellar rotation, the peak disk 



accretion rate is su fficiently high to ensure effec tive neutrino 
cooling of the disk l|Chen fc Beloborodovll2007l ). 



3.3 Polytrope model 

Finally, we consider the model polytrope with index n = 3, 
which could be used to describe the cores of most mas- 
sive s tars at the pre-supernova phase l|Tutukov fc Fedoroval 
2007). In this model the concentration of mass towards the 
centre is much weaker, resulting in higher moment of iner- 
tia and larger angular momentum compared to the Bethe's 
model with the same mass, radius and rotation frequency. 
Even if we consider models with the same spin parameter 
a s , the polytrope yields generally higher fraction of mass 
accreted via the accretion disk and more rapidly rotating 
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Figure 6. The evolution of black hole's mass and spin in the polytrope model of progenitor with the initial black hole mass M c = M s /9. 
The left panel shows the total mass of the black hole (thick solid line) as well as the mass accumulated via the accretion disk for different 
rotation rates of the progenitor, a s = J a c/GM 3 = 0.33, 0.58, 1.0, 1.7, 3.0 (dash-dotted lines). Faster rising lines correspond to higher 
rotation rate. The right panel shows the spin parameter of the black hole, a, for the same values of a s ■ 
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Figure 7. The evolution of black hole spin for polytropic models 
with a s = 3 and M C /M B = 1/3 (dot-dashed line), 1/5 (dotted 
line), 1/9 (dashed line), and 1/31 (solid line). 

black holes (see Figure [6]). Similarly to other models, the 
final value of the black hole's spin does not show strong 
dependence on the iron core mass fraction, at least for 
M c /M s € (1/3, 1/31) (see Figure^. 

The polytrope model was also used to test our calcu- 
lations against the fully general relativistic simulations by 
IShibata fc Shapiro] 1|2002| ). For the polytropic star with an- 
gular momentum a 3 — 1 our model gives a black hole with 
Mi, = 0.90M S and a = 0.76 by the time of disk formation. 
This is in excellent agreement with the numerical simula- 
tions which give M& = 0.90 and a = 0.75. 



4 JET SIMULATIONS 

The analysis carried out in Sections and [3] suggests that 
during the collapse of a WR star in a very close binary sys- 
tem, the conditions can become favourable to production of 



LGRB jets either via the neutrino heating or the Blandford- 
Znajek mechanism. Although the production of jets via the 
Blandford-Znajek mechanism has already been studied nu- 
merically in several previous papers the conditions suggested 
by the binary scenario are different from those explored so 
far. By the time of the accretion disk formation the black 
hole is much more massive compared to the usually assumed 
Mh — 2Mq. Its rotation rate is noticeably lower compared 
to a ~ 1, assumed in the past. Finally, the progenitor's rota- 
tion is not differential but uniform. These differences invite 
additional numerical simulations to explore the new region 
of parameter space. 

4.1 Setup of Simulations 

The progenitor model describes a compact WR star of radius 
R s = 3 x 10 10 cm and rotation period T s = 1.4hr; the cor- 
responding specific angular momentum on the stellar equa- 
tor is l s — 1.13 x 10 18 cm 2 s _1 . The progenitor's magnetic 
field is assumed to be purely poloidal and uniform, with the 
strength B = 1.4 - 8.4 x 10 7 G. 

Simulations of this type are computationally expensive 
even in 2D. On the other hand, the early stages of the col- 
lapse are very simple and can be treated analytically with 
sufficient accuracy. For these reasons, we start simulations 
only after the expected time of the disk formation, t 3 = 17s. 
Based on the analysis given in the previous sections, the 
black hole mass is set to Mh = IQMq and the mass accre- 
tion rate to 0.14Mqs _1 . The initial radial distributions of 
mass and velocity are the same as in the Bethe model: 

p oc R" 3/2 , v r = {2GM h /R) 1/2 . (32) 

The initial distributions of angular momentum and magnetic 
field are derived from the progenitor distributions by taking 
into account the distortions caused by the free-fall collapse 
over the time t s : 

( t \ 4/3 

l(R,6) = n 3 (Rsm6) 2 1 + — , (33) 
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Table 1. Numerical models. 



Model 


a 


Bo 


*28 


Lbz 


Mh 


L B z/M h c 2 


Ml 


0.6 


1.4 


0.46 








M2 


0.6 


4.2 


1.5 


0.44 


0.017 


0.0144 


M3 


0.45 


8.4 


3.1 


1.1 


0.012 


0.049 



a is the black hole spin; Bo is the initial magnetic field strength 
in the units of 10 7 G; \E"28 is the magnetic flux accumulated by 
the black hole by the time of explosion in the units of 10 28 G cm 2 ; 
Lbz is the total power of the BZ-mechanism during the explosion 
in the units of 10 51 erg/s; is the black hole mass accretion 
rate during the explosion in the units of Mq s . 



r Bo sin 6 cos ( 
B = 



V7 



R 2 (l + t/t ff (R)f 3 : 



B 



„ Bo sin 2 6 
= 2^7 



2R{l+t/tf f (R)) 1/a : 



(34) 
(35) 



where t ff (R) = y // 2R 3 /9GMh is the local free fall time scale 
and 7 is the determinant of the metric tensor of space (see 
Appendix[X]| . In these simulations we studied three different 
cases summarised in Table [T] 

The simulations were c arried out with 2D axisymm etric 
GRMHD code described in iKomissarovl |l999, 2004). The 
gravity effects are introduced via the Kerr metric with fixed 
parameters; the Kerr-Schild coordinates are used in order to 
avoid the coordinate singularity at the horizon. The compu- 
tational grid is uniform in the polar angle, 9, where it has 
180 cells and logarithmic in the spherical radius, R, where 
it has 445 cells. The inner boundary is located just inside 
the event horizon and adopts the free-flow boundary condi- 
tions. The outer boundary is located at R — 8.3 x 10 9 cm 
and at this boundary the flow is prescribed according to the 
Bethe's model. 

In the simulations we used realistic equation of state 
(EOS) that takes into account the contributions from radia- 
tion, lepton gas (including pair plasma) , and non-degenerate 
nuclei (hydrogen, helium, and oxygen). This is achieved via 
incor poration of the EOS code HELM l|Timmes fc Swestvl 
2000), The neutrino cooling is computed assuming opti- 
cally thin regime an d takes into account URCA-processes 
|lvanova et al.l Il969h . pair annihilation, p hoto-production, 
and plasma emission (ISchinder et al. 1987) , as well as syn- 
chrotron neutrino emission (|Bezchastnov et aL | Il997t ). In 
fact, URCA-processes strongly dominate over other mech- 
anisms in this problem. Photo-disintegration of nuclei is in- 
cluded via m odification of t h e EO S following the prescrip- 
tion given in lArdelian et all |2005). The equation for mass 
fractio n of free nucleons is adopted from lWooslev fc Baronl 
(|1992| ). We have not included the radiative heating due to 
annihilation of neutrinos and antineutrinos produced in the 
accretion disk mainly because this requires elaborate and 
time consuming calculations of neutrino transport. 



4.2 Results 

In general, the results of these si mulations are in agree- 
ment w ith our previous studies ([Barkov fc Komissarov I 
l2008al lbl; lKomissarov fc Barkov \\200$ ). Because of the mod- 
ified setup, which corresponds to the later stages of the col- 
lapse, the accretion disk if formed straight away. At the same 



time, the accretion shock separates from the disk surface and 
quickly expands up to R ~ 100 — 200i? 9 . In the model Ml, 
the shock then begins to oscillate and no jets emerge by the 
end of the simulations, t — 19.5s. In contrast, both the mod- 
els M2 and M3 eventually develop polar jets of relativistic 
plasma which are powered via the Blandford-Znajek mech- 
anism (see Figs(8]and|9]). These results comply with the BZ 
activation condition, k > 0.2, where 



47rr £ 



(36) 



is the activation parameter, $ is the magnetic flux thread- 
ing the black h ole and M is the mass accreti on rate of the 
collapsing star (Komiss arov fc Barkov ||2009T ). For the pa- 
rameters of the present simulations we have k ~ 0.07^28, 
where the magnetic flux is given in the units of 10 28 Gcm 2 , 
and, thus, one would expect the BZ mechanism to become 
activated for '3/28 > 3. As one can see from the data pre- 
sented in Table 1 which was indeed the case. 

According to the simple monopole model of black hole 
magnetosphere, the power of the BZ mechanism is 



S BZ = 1.4xl0 6 V 2 (a)^ 8 ( I ^- 



ergs 



(37) 



where /2(a) = a 2 (l + y/l — a 2 ) 2 (jBarkov fc Komissarov I 
2008a). Like in our previous simulations, the direct measure- 
ments of energy flux across the BH horizon roughly agree 
with this result (see Table 1). 

One significant difference with the results of previ- 
ous simulations is the development of a one-sided jet in 
model M2. Although noticeable deviations from the equato- 
rial symmetry have been observed before, in particular the 
asymmetric oscillations of the accretion shock, such a strong 
deviation is observed for the first time. The initial solution is 
not exactly symmetric because of the rounding errors, but 
they are tiny and the observed braking of the equatorial 
symmetry has to be rooted in the nonlinear dynamics of the 
flow. It appears that the accretion flow, which is deflected to- 
wards the equatorial plane at the oblique shock driven by the 
northern jet, protrudes into the southern hemisphere. There 
it collides with the accretion flow of the southern hemisphere 
and together they stream towards the black hole's southern 
pole, thus suppressing the development of a southern jet. 

If persistent, such a one-side jet could impart a strong 
kick on the black hole and the bin ary, significantly a ltering 
its motion in the parent galaxy (|Fragos et al. 2009). The 
maximum kick velocity can be estimated as 



E 



jet 



cM h 



170 



E 



jet 



10 52 erg I V M h 



(38) 



which is consistent wi th the observations of the X-ray bi- 
nary XTEJ1118+480 (|Gualandris et al.ll2005l ). However, at 
present we cannot say whether such one-sidedness can per- 
sist during the life-time of LGRB or this is just a transient 
phenomenon. 



5 MERGER SCENARIO 

The case of close tidally-locked binary considered above in- 
volves binaries with orbital separation very close to the size 
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Figure 8. Model M2 (Bo = 2.2 X 10 7 G, a = 0.6) at the time of Figure 9. Model M3 (Bo = 8.8 X 10 7 G, a = 0.45) at t = 1.35s 

t = 1.35s after the start of simulations (18.3s after the start of the after the start Q f simulations (18.3s after the start of the collapse), 

stellar collapse). The colour image shows the baryonic rest mass The colour image shows thf , baryomc rest mass densityi l O g 10 p, 

density, logiop, in CGS units, the contours show the magnetic in CGS unitSj thc cont0 urs show the magnetic field lines, and the 

field lines, and the arrows show the velocity field. arrows show the velocity field 



of the Roche lobe of the WR star and this suggests to go one 
step further and consider the case of even smaller separation 
which can lead to the com mon envelope evolution resu lting 
in a merger of t he binary dTutukov fc Yungelsonll 19791 ) and 
GRB explosion (|Zhang fc Frverll200l[ h Such a merger can be 
divided into tree phases. During the fist phase, the compact 
companion spirals inside the extended envelope of the nor- 
mal star and spins it up via deposition of its orbital angular 
momentum. The compact star also increases its mass and 
spin via th e Bondi-type a c cretio n. According to the sim- 
ulations of IZhang fc Frvej (120011 ) , during the last 500s of 
the in-spiral the compact star can accumulate up to 3.5Mq. 
Thus, the mean accretion rate is less than 10 _2 Mq/s imply- 
ing inefficient neutrino heating. 

The second stage begins when the compact star ap- 
proaches the c entre of its WR comp anion and the accretion 
rate increases. IZhang fc Frve3 (|200lf) find that in the case of 
I6M0 companion the neutrino annihilation mechanism can 
operate for around 60 seconds and release about 10 52 erg. 
This is more than enough to drive a supernova explosion. 
For the companion mass below 8Mq the neutrino heating is 
too weak and the second phase is absent. 

The third phase takes place if the second phase does 
not result in the supernova explosion or if the explosion is 
highly non-spherical and does not remove the equatorial lay- 



ers of the WR star. During this phase the compact object, 
already a black hole, accretes these layers, which have been 
spun up during the first phase. Assuming that the mass of 
the compact star is small compared to the mass of its WR 
companion, its orbital angular momentum can be found via 
the Keplerian law 

J C (R) = M c yJ GM(R)R, 

where M(R) is the WR mass inside the radius R. As the 
compact star moves from the radius R to R + dR, it trans- 
fers the angular momentum dJ c (R) — (dJ c /dR)dR to the 
envelope of the WR star. Assuming that most of this angu- 
lar momentum is transferred to the mass dM = (dM / dR)dR 
of the envelope located between R and R+dR, we obtain the 
specific angular momentum of the envelope after the merger 
as 

dj c _ dJc/dR 
~~dM~ dM/dR' 

For the Bethe's model, where M(R) is given by Eq. (|18[) . this 
gives 

,K x (jm) 1/2 ( 1+HR/R °»' (39) 

which is smaller than the local Keplerian angular momentum 
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provided M(R) > M c (l + \nR/R c )/2. This suggests that if 
M B ^> M c then only a small fraction of the common envelop 
is lost during the merger. For R — R s this equation gives 



I ~ 5.2 x 10 



M c 

2Mr. 



R s 

Rq 



1/2 



( Ms 



\WM t 



-1/2 



-\(40) 



In the a-model, the accretion time scale of the disk with 
such angular momentum can be estimated via 

td 



: 8000 s 



ao 
OOl 



a<5 2 (GMs) 2 

3/2 



Rs_ 

Rq 



M 3 



\2MqJ \10Mq 



-7/2 



(41) 



This is significantly longer than the duration of the stellar 
collapse (see Eq |15p . In fact, such a long time scale suggests 
the possibility of explaining the phase of shallow decay and 
late fl ares in the X - ray light curves of LG RBs discovered by 
Swift (jZhang||2007l ; IChincarini et alj|2007h . 

To find the mass accretion rate as a function of time we 
note that 



M = 



dM _ dM/dR 
did did I dR 



Using Eqs. (|18l41l39|) to evaluate dM/dR and dt d /dR we ob- 
tairQ 



3ln(R 3 /Rc) t 



Thus, on the time scale of 10 — 10 s the mass accretion rate 
is very low, M ~ 1O~ 3 4TO _4 M0S _1 , ruling out the neutrino 
mechanism and leaving the BZ mechanism clear favourite. 
Indeed, the maximum possible amount of magnetic flux that 
can be accumulated by the black hole is given by the balance 
of magnetic pressure and the gas pressure of the accretion 
disk, 



8 Ti 



PC a 



(43) 



where c a is the sound speed. If we utilize the model of a-disk 
and estimate the magnetic field strength at the gravitational 
radius then the corresponding magnetic flux will be 



3 x 10" 



M b \ 



(^-Y 1/2 (^^ iM^Gcm 2 , (44) 
VO.03/ \WM Q J 1 



where M\ is the mass accretion rate in the units of Mq/s. 
Even for Mi as small as 10 -4 this equation gives the sub- 
stantial value of $„, ~ 3 x 10 27 Gcm 2 . The corresponding 
BZ power, E BZ ~ 2.2 x 10 49 erg/s, is more than sufficient to 
explain the X-ray observations, allowing the magnetic field 
to be even weaker compared to the value suggested by Eq|43l 



6 DISCUSSION AND CONCLUSIONS 

One of the main issues of this study was to investigate the 
efficiency of the tidal spin up in close massive binaries in the 
context of the collapsar model of LGRBs. In particular, we 



7 This equation is the same as Eq |31l but t spans a different range 
of timescales, now dictated by the disc accretion time. 



wanted to find out the typical masses and spins of the black 
holes formed during the collapse of WR companion. It turns 
out, that the BH spin in this model is rather modest. For 
example, in the most optimistic case of a binary with the 
smallest possible orbital separation, the spin parameter of 
the WR star is relatively high, a a — 6, and one may have 
expected the BH to be rapidly rotating. However, we find 
that the spin parameter is only a ~ 0.4 at the time of the 
accretion disk formation, and a ~ 0.8 by the end of the stel- 
lar collapse, which is significantly lower than the maximally 
possible value a = 1. This is mainly due to the significant 
loss of angular momentum suffered by the mass accreted via 
the disc; the rejected angular momentum is either stored 
in the remote part of the accretion disk or removed by a 
disk wind. Indeed, as soon as the accretion disc is formed, 
the rate of accretion of angular momentum slows down sig- 
nificantly. Moreover, by the time of the disk formation the 
black hole mass is already rather high, exceeding half of the 
progenitor mass prior to the collapse. Thus, the black hole 
simply runs out of accreting matter bef ore its rotatio n can 
approach the maximal possible rate (~cf. lThorneHl974T ). 

The mass accretion rate in this scenario is much lower 
compared to the usual M = (0. 1 4- 1)Mq s -1 invoked 
in th e standard collapsar model (|MacFadven fc Wooslevl 
1999). This makes the neutrino mechanism less attrac- 
tive compared to the magnetic mechanisms, and the BZ- 
mechanism in particular. In fact, the very rapid decline in 
the efficienc y of neutrino mecha n ism below M < 0.02 



0.05A fo s' ^Popham et al.1 Il999l ; IZalamea fc Beloborodovl 



I 20091 ) makes the explanation of the LGRB bursts with du- 
ration > 100s rather problematic even within the standard 
collapsar model due to the low mass accretion rate expected 
on such time scale (see Eg 1421 and Fig[5}. 

However, the BZ mechanism could have its own difficul- 
ties in this scenario. Indeed, it requires very strong ordered 
magnetic field. For example, in order to provide the power of 
lO^ergs -1 the black hole of mass WMq and a = 0.6 should 
accumulate the magnetic flux of order f = 8 x 10 27 Gcm 2 . 
The magnetic flux necessarily to activate the BZ mechanism 
soon after the formation of accretion disk is even higher. 
According to Table I this is of the order fewx 10 28 G cm 2 . 
Perhaps the free-fall accretion rate set up in our simulation 
is a bit too high and could have been reduced by a factor 
of 10. However according to Eq. (|36[) . this would reduce the 
critical value of magnetic flux only by a factor of 3. 

The origin of such strong field is not clear. It 
could be gene rated via magnetic dyn amo in the accre- 
tion disc (e.g. [Brandenb urg et al.l I1995T) or in the convec- 
tive c ore of the progenitor (e.g. ICharbonneau fc~M acGrcgor 
2001). It may also be inherited by the progenitor from 
the interstellar medium (I SM) during its formation (e.g. 
iBraithwaite fc Spruit! 120041 '). The current status of both 
the stellar and disk dynamo theories does not really al- 
low to make reliable conclusions. Even the issue of ad- 
vection of externally generated magnetic field by the ac- 
creti on disk onto the central black hole is still unresolved 
(e.g. Bisnovatyi-Kogan & Ruzmaikin 1974; van Ballcgooiien 
' 19891; iGuan fc Gammiel 120091; I Spruit fc Uzdenskyj 120051 ; 
Rothstein fc Lovelace! 120081 ; llgumenshchev I l2008h . There 



seems to be a general agreement that accretion disks pro- 
duce mainly azimuthal magnetic field and unable to gen- 
erate poloidal field on scales exceeding the disk hight. The 
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magnetic dynamo in convective cores of massive stars could 
be more promising in this respect. For exa mple, from the 
results of lCharbonneau fc MacGreeorl |2001) it seems possi- 
ble to generate up to $ ~ 10 28 Gcm 2 in the convective cores 
of B stars. 

By design, our numerical simulations cannot address 
the issue of magnetic field generation is accretion disks and 
strictly-speaking deal only wit h the fossil model of magnetic 
field. The numerical results by iBraithwaite fc Spruitl (I2004T ) 
suggest that strong fossil field can relax to a simple ordered 
configuration with dipolar poloidal field on a relatively short 
time scale, which makes our setup not that unrealistic. How- 
ever, further studies are required to verify this model. The 
observations of massive stars do not support magnetic flux 
of order 10 28 G cm 2 and higher. The current record is held by 
1 Ori C, whose dip olar magnetic flux * ~ 2 x 10 27 Gcm 2 , 
Donati et all l|2002f n. One may speculate that most of the 
magnetic flux is hidden in the stellar interior. Indeed, the 
resistive time scale across the extended radiative outer lay- 
ers of massiv e stars exceeds their life tim e by many orders 
of magnitude IBraithwaite fc Spruitl (120041 ). 

The fact that the magnetic flux of neutron stars is less 
than 10 27 Gcm 2 also seems to be working against the fossil 
hypothesis. However, neutron stars are collapsed compact 
Fe cores of massive stars. The typical cross section of such a 
core is several orders of magnitude below that of the whole 
star and, thus, the core may account only for a small frac- 
tion of the total magnetic flux hidden inside the supernova 
progenitor. 

The host galaxies of LGRBs show strong evidence of 
enhanced star formation ( Bloom et al.l 1 19981 ; I Sokolov et al.l 
l200ll ; iFruchter et al ] |2006f T It is interesting that the recent 
observations of such starburst galaxies also indicate strong 
ISM magnetic field, in fact up to ten times stronger com - 
pared to the Milky Way [|Beck fc Krausdl 20051 : iBeckll 20081 ). 
This suggests that magnetization of young stars in the host 
galaxies of LGRBs can be abnormally high as well. 

Another interesting proposal ste ms from th e theor y of 
Sun's magnetic activity proposed by lUzdenskvl l|2007l ). In 
particularly, he argued that fast reconnection can only op- 
erate in collisionless plasma and in the collisional regime 
the reconnection rate reduces to the much slower rate of 
Sweet-Parker. Since the collapsar plasma is collisional even 
in the rarefied funnel of the accretion disk then, accord- 
ing to this theory, the reconnection rate in the black hole 
magnetosphere can be relatively slow. An additional unex- 
plored factor in the LGRB context is the effects of quan- 
tum physics on magnetic reconnection. Indeed, the expected 
magnetic field strength is well above the quantum value of 
B q = m 2 c 3 /fte = 4 x 10 13 G. One may speculate that under 
this conditions the reconnection rate becomes even slower. 

In the case of slow reconnection, the black hole may be 
able to build strong magnetic field via collecting the alter- 
nating magnetic field generated in the accretion disk. Since 
the magnetic stresses are invariant with respect to change 
of magnetic polarity such striped structure of magnetic field 
has no effect on the efficiency of the BZ-mechanism. Fur- 
ther downstream of the LGRB flow, where its plasma be- 
comes collisionless or the magnetic field becomes sufficiently 
weak, the reconnection accelerates. However, as long as this 
occurs beyond the Alfven surface, which for a black hole 
with reasonable spin does not greatly exceed the gravi- 



tational radius, this does not disrupt the near magneto- 
sphere of the black hole and does not reduce the efficiency 
of the BZ-mechanism. Moreover, such delayed reconnec- 
tion could promote bulk acceleration of the LGRB flow 
|Drenkhahn fc Spruitlliooll . 

Finally, the neutrino heating of the polar region, not in- 
cluded in our analysis and simulations, may also play a very 
important role, by initiating the LGRB outflow and creat- 
ing the low density chanel in the polar direction early on, 
when the mass accretion rate is still sufficiently high for ef- 
fective neutrino-antineutrino annihilation. This would allow 
the BZ-mechanism to be activated along the field lines fill- 
ing the channel even if the black hole magnetic flux is much 
lower compared to the values quoted above. Later on, when 
the mass accretion rate drops and the neutrino mechanism 
can no longer provide sufficient power, the BZ-mechanism 
can take over the role of main driver of the LGRB flow. 
One may even contemplate the scenario where the GRB pre- 
cursors are related to the neutrino-driven stellar explosions 
and the main bursts to the magnetically-driven BH jets un- 
leashed in the space cleared up by the blast. The delay be- 
tween the two phases could be related t o the disruption of 
the a ccretion flow by the supernova blast (jWang fc Meszaroa 
120071 ). Because of the rotational effects the disruption may 
not be as severe in the equatorial direction, compared to the 
polar direction, a s in t he one-dimensional simulations by 
iMacFadven et all (|200ll ). leading to shorter fallback time- 
scales. The magnetic jets, though very powerful, could be 
less disruptive compared to the neutrino-driven jets because 
the magnetic hoop stress, associated with the azimuthal 
component of magnetic field, makes the sideways expansion 
of the jet cocoon less effective. 

The most interesting, in view of the recent Swift ob- 
servations of LGRB afterglows, version of the close binary 
scenario for GRB progenitors is the common envelope case, 
where the compact star, either a black hole from the on- 
set or a neutron star which eventually collapses into a black 
hole, spirals inside the normal WR star. The large angular 
momentum transferred to the external layers of the WR star 
quite naturally leads to long accretion time scales, ~ 10 4 s. 
Thus, the central engine of LGRB jets arising in this sce- 
nario could operate for a sufficiently long time to explains 
the sh allow phase of the X-ray light curves discovered by 
Swift (|Zhanell2007l ). The X-ray flares, which are often seen 
during this phase, may result from the gravitat ional insta- 
bilities developing in this disc (jPerna et al.l l2006). Although 
the BZ mechanism is not that sensitive to the mass accretion 
rate as the neutrino mechanism, some dependence is still 
expected. For example, Eqs. (l37l44p suggest that the power 
of the BZ mechanism may be proportional to the mass ac- 
cretion rate. This can explain why the gamma ray emission 
becomes undetectable on the time-scale of the shallow decay 
of X-ray afterglows. 

The extremely high rotation rates, about 50% of the 
br eak-up speed, a s sumed in the single progen i tor m odel 
bv lYoon fc Langerl (|2005l ) and lWooslev fc Hegeri (|2006l ) im- 
ply that in this model the outer layers of the collapsing 
star can also develop long-liv ed accretion disk. Indeed , in 
the "showcase" model 16TI of lWooslev fc Hegeri l|2006l ) the 
outer ~ 2Mq have the specific angular momentum increas- 
ing outwards from 10 18 to 10 19 cm 2 s -1 at the pre-supernova 
phase. According to Eq |40l this corresponds to the disc ac- 
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cretion time scales of order 10 4 s. However, such a long time 
scale still rules out the neutrino annihilation as the mecha- 
nism for powering the collapsar jets. 
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l; , _ Bp sin 6 cos 9 ^ . ^ 



V7 



and 



B = 



B sin 2 . 



V7 



tff(R) 



1/3 



4/3 



R 1 + 



tff(R) 



(Al) 



(A2) 



where 7 is the determinant of the metric tensor of space 
and the vector components are given in the non-normalised 
coordinate basis, d/dx 1 . This approa ch has been used in 
Bisnovatyi-Kogan & RuzmaikinI (| 19741 ). 



APPENDIX A: EVOLUTION OF ANGULAR 
MOMENTUM AND MAGNETIC FIELD IN THE 
BETHE'S MODEL OF STELLAR COLLAPSE 

The free fall model bv lBethel (|l990h approximates the kine- 
matics of stellar collapse by the model 

dR _ ( if t < 0, 

dt ~ \ -(2GM(R)/R)- 1/2 if t > 0, 

where R is the radius of collapsing shell and M(R) is the 
mass inside this radius. Since dM/dt = 0, this equation is 
easily integrated 

Ro(R,t) = R(l + t/t ff (R)) 2/3 = const, 

where R = R(0) and t ff (R) = V / 2R 3 /9GM(R) is the local 
free fall time. 

Given the initial distribution of angular momentum, 
lo = Q.(Ros'm6) 2 , the conservation of angular momentum 
yields 

l(R,t) = l (Ro(R,T)) = fi s (i?sin#) 2 (l + t/t ff (R)) i/3 . 
Similarly, the conservation of magnetic flux requires 

V(R,e,t) = V (R {R,t),6). 
For the uniform initial magnetic field, 

Vo(Ra,0) = vrBosin 2 ^- 

Thus, 

$(R,6,t) = TTB sm 2 6R 2 (l + t/t ff (R)) 4/3 . 
The poloidal magnetic field can be found via 

B i = — e* jV 9,-* 
p 2tt 3 ' 

where e ljfc is the Levi-Civita tensor of space. This gives us 
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